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ABSTRACT 

The lack of detected pulsars at the Galactic Center (GC) region is a long-standing mystery. We 
argue that the high stellar density in the central parsec around the GC is likely to result in a pulsar 
population dominated by millisecond pulsars (MSPs), similar to the situation in globular cluster 
environments. Earlier GC pulsar searches have been largely insensitive to such an MSP population, 
accounting for the lack of pulsar detections. We estimate the best search frequency for such an 
MSP population with present and upcoming broad-band radio telescopes for two possible scattering 
scenarios, the “weak-scattering” case suggested by the recent detection of a magnetar close to the 
GC, and the “strong-scattering” case, with the scattering screen located close to the GC. The optimal 
search frequencies are ~ 8 GHz (weak-scattering) and ss 25 GHz (strong-scattering), for pulsars with 
periods 1 — 20 ms, assuming that GC pulsars have a luminosity distribution similar to that those in 
the rest of the Milky Way. We find that 10 — 30 hour integrations with the Very Large Array and the 
Green Bank Telescope would be sufficient to detect MSPs at the GC distance in the weak-scattering 
case. However, if the strong-scattering case is indeed applicable to the GC, observations with the full 
Square Kilometre Array would be needed to detect the putative MSP population. 

Subject headings: Galaxy: centre - pulsars: 


1. INTRODUCTION 

The allure of testing General Relativity near a super- 
massive black hole and in vestigating the latter’s accre¬ 
t ion e nvironment (see, e.g JPfahl fc Loebll2004 [Liu et all 
120121) has provided the motivation for numerous searches 
for radio pulsars at the Galactic Center (GC) over the 
last two decades (e.g. Johnston et al.ll2006t iDeneval 120101 
iMacauart et alJ I201CIL Yet, despite predictions that 
100 —1000 radio pul sars with orbital periods < 100 years 
are o r biting Sgr A* dPfahl fc Loebl 12004 iWharton et al. 
2012; IChennamangalam fc Lorimeil 12014 IZhang et al. 
2014|), not a single normal pulsar has so far been detected 
within 10' of Sgr A*. 

A long-standing problem related to these pulsar sur¬ 
veys relates to the scattering environment at the GC. 
The extreme angular broadening exhibi ted by Sgr A* 
and OH/IR star s in the GC region (e.g. ILo et al.l 119851 
iFrail et all 1 1 9911) . coupled with inferences about the dis¬ 
tribution of the scattering material, led to the conclusion 
that the pulse smearing caused by turbulence along the 
line of sight to the GC is extreme, with a temporal smear¬ 
ing time of as 10 1 2 3 v^„ s, where is the observing 

frequ ency in GHz (jCordes fc Laziolfl997t lLazio fc Cordesl 
Il998f) . Such a large temporal smearing time would ob¬ 
struct the detection of pulsed radio emission at the usual 
low frequencies at which pulsar searches are usually car¬ 
ried out (due to their steep radio spectrum). Indeed, 
for normal pulsars, with periods ss 0.5 s, the tempo¬ 
ral smearing w ould become negligible only at frequencies 
above 15 GHz llMacauart et al.ll2010l h where the steep¬ 
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ness of the pulsar emission renders them weak and hence 
difficult to detect, even with today’s most sensitive tele¬ 
scopes. 

The recent detection of the GC magnetar PSRJ1745- 
29, located within 3" of the p osition of Sgr A* 
(|Kennea et al.ll2013i: IMori et al.ll2013l h has caused an up¬ 
heaval in the field, lending ren ewed im petus t o sear ches 
for pulsars close to Sgr A*. iSpitler et alJ (120141) re¬ 
ported the detection of pulsed emission from the magne¬ 
tar at frequencies as low as 1.1 GHz, inferring a temporal 
smearing timescale of only (1.3 ± 0.2)^qj^ z s. The effect 
of pulse smearing on at least one object in the GC envi¬ 
ronment is thus three orders of magnitude smaller than 
that inferred from earlier studies. The magnetar also ex- 
hi bits the same deg ree of angular broadening as Sgr A* 
(jBower et al.ll20l4 ). suggesting that it lies behind the 
same hyperturbulent scattering region. This hypothesis 
is sup ported by the high rotation measure of the mag - 
netar (|Eatough et al.lDOTjl : iShannon fc JohnstonlfeOlfili . 
which indicates that its proximity to Sgr A* on the sky 
is not merely a chance alignment. 

The detection of the magnetar close to the GC sug¬ 
gests that the pulse smearing towards the GC is also 
relatively benign, assuming that the scattering screen 
is uniform. This would imply that radiation from or¬ 
dinary pulsars at the GC should be readily detectable 
at frequencies > 3 GHz. However, if the effects of 
pulse smearing are similarly small over the entire GC 
region, previous surveys should have detected a signif¬ 
icant fraction of the pulsar d i stribu t ion with “normal” 
spin periods (Wharto n et al.l 120121 IDexter fc O’LearvI 
l20H . This dearth of pulsar detections is particularly 
acute given that magnetars ar e believed to represent 
only ~ 0.2% of all radio pulsars (lOlausen fe Kaspill20l4 

IChennamangalam fe Lorimeil 120141) . The lack of “slow” 
pulsar detections suggests that they either constitute an 
unusually small fraction of the entire GC pulsar popula- 
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tion, or are significantly under-luminous relative to the 
population of slow pulsars that have been detected else¬ 
where in the Galaxy. 

In p assing, we note that lChennamangalam fe Lorimed 
( 2014 1 used a Bayesian analysis combined with an as¬ 
sumed log-normal pulsar luminosity function to find that 
existing pulsar surveys at the GC are not sufficiently 
deep to eliminate the possibility that a substantial popu¬ 
lation of non-recycled low-luminosity pulsars might exist 
at the GC. While this is broadly correct, the shape of the 
pulsar luminosity function is not at all well-constrained 
at the low-luminosity end. As such, there is consider¬ 
able error in the extrapolation of the distribution to low 
luminosities, and hence on such estimates of the total 
size of the normal pulsar population in the central par¬ 
sec. For example, an uncertainly of only 20% in t he two 
parameters of the log-normal distribution (jBagchi et, al.1 
120111 : iChennamangalam fc Lorimed 12014 (0 yields an un¬ 
certainty of two orders of magnitude in the survey com¬ 
pleteness. Similar estimates of the number of potential 
MSPs at the GC are plagued with an even greater degree 
of uncertainty because the luminosity function of recy¬ 
cled pulsars in globular clusters is even less well known. 
In the present analysis, we will adopt a pragmatic ap¬ 
proach to the pulsar luminosity distribution of pulsars: 
our detection arguments are made only with reference to 
the known luminosities of detected Galactic pulsars. 

In this paper, we explore the reasons for the paucity 
of pulsar detections at the Galactic Center. In §E1 we 
discuss the nature of the pulsar population at the GC re¬ 
gion, and argue that it is likely to be dominated by recy¬ 
cled millisecond pulsars (MSPs); present pulsar surveys 
of the GC have been largely insensitive to such a pop¬ 
ulation. Next, we present in §[3] detailed calculations of 
the expected signal-to-noise ratio for a wide-bandwidth 
search for MSPs at the Galactic Center with present and 
future telescopes, incorporating the frequency depen¬ 
dence of temporal smearing and sky temperature across 
the different observing bands. Finally, the results of this 
work are summarized in §01 


2. A MILLISECOND PULSAR POPULATION AT THE 
GALACTIC CENTER 

Over the last two decades, several searches have been 
carried out for pulsars at the GC, mostl y at frequencies 
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higher observing frequencies are more amenable to the 
discovery of pulsars, with the optimal frequency range 
for searches for “normal” pulsars - those with periods 
of ss 0.5 seconds being 10 — 16 GHz. Following this, 
there have been a number of deep searches at f requencies 
above 10 GHz, at 15 GHz by 
12 — 18 GHz bv ISiemion et al 


Macquart et al. (1201011 , at 
(1201311 . and at 19 GHz by 
lEatough et al.1 (12013H . Remarkably, despite integration 
times exceeding 10 hours with 100-m single-dish tele¬ 
scopes, none of these searches has discovered a single 
pulsar in the Galactic Centre region! 

Using the tempo ra l smearing e s timate s of 
lLazio fc Cordesl (119981) . iMacouart et al.1 (1201011 esti¬ 
mated that their 15 GHz search would have been 
sensitive to ~ 15% of the Galactic center pulsar 
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Fig. 1.— The 1.4 GHz luminosity (in m jy kpc 2 ) of the known pul¬ 
sar population lIManchester et al.120051 ) is pl otted versus pulsar pe - 
riod. The 1 0<x sensi tivities of previous 5 GHz (Johnston ct a l .120061 ) 
and 15 GHz GBT UMacauart et al1l2010l) searches of the Galactic 
Center are shown by the green and red curves, respectively, with 
the dashed and solid curves representing the “strong” and “weak” 
temporal smearing scenarios, respectively. 

population, assuming the luminosity distribution of 
known Galactic pulsars. However, if the pulse smearing 
is benign, as suggested by the detection of the GC 
magnetar, then even the earlier lower-frequency searches 
would have been sensitive to normal pulsa rs located near 
the GC. Indeed. [Dexter fc O’LearvI (120141 ) estim ate that 
both the 5 GHz sea rch ofIJohnston et al.1 ( 2006fl and the 
15 GHz search of iMacauart et al.l (I201C ) would both 
have been sensitive to ~ 20% of the GC pulsar popula¬ 
tion, if its characteristics resemble t hose of the known 
pulsa r population (see Fig. 2 of iDexter fc O’LearvI 
Soil. This “mi ssing pulsar” problem was used by 
Dexter fc O’LearvI ([20141 ) to argue that the GC pulsar 
population may be dominated by magnetars, i.e. that 
the population is very different from that in the rest 
of the Galaxy (where, as noted earlier, only ss 0.2% of 
known radio pulsars are magnetars, although the mag¬ 
netar birth rate may be 10 — 50% of the total neutron 
star birth rate in the Galaxy; e.g. iKeane fc Krameil 

HHI). 

However, a number of studies have raised another pos¬ 
sibility, namely that the GC pulsar population is domi¬ 
nated by recycled millisecond pulsars (MSPs). The dense 
stellar environment at the GC is likely to result in spin¬ 
ning pulsars up to millisecond periods by frequent close 
interactions with neighbouring stars, by analogy with 
the population of MSPs detected in the dense ste llar en¬ 
vironments of globular clusters ( e .g. lAlpar et al. Il982t 
iVerbuntl 119871 : iCamilo et alJ 120001 : iRansom et al.N2005D . 
Note that the stellar density in the central parsec of 
the GC is ~ 10 6 per cu bic parsec (e.g. iGenzel et alJ 
119961 : iSchodel et al i [20071) . a couple of orders of mag¬ 
nitude larger than the stellar density in globular cluster 
cores (< 10 4 per cubic parsec), implying that close in¬ 
teractions are far more likely in the vicinity of the GC. 
The formation rate of low-mass X-ray binaries (LMXBs) 
in dense environments is also proportional to both the 
num ber density of n eutro n stars and the stellar den¬ 
sity (lYerbunt fc Hudfl987l) . Since the neutron star den¬ 
sity itself scales with the stellar density, this implies 
that the formation rate of LMXBs (and hence, that of 
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binar y pulsars; e.g. lAlpar et~ahl 119821 : iCampana et al.1 
1998) is roughly proportional to the square of the stel¬ 
lar density. When coupled with the higher stellar den¬ 
sity at the GC, this implies a far higher formation rate 
of LMXBs (and binary pulsars) than in globular clus¬ 
ter cores. We note that more than 90% of all known 
globular cluster pulsars have periods smaller than 30 ms 
((Manchester et al] 1200511 . indicating that a dense envi¬ 
ronment can dramatically alter the period distribution. 
Further, an over-abundance of X-ray transients has been 
detected within 1 pc of Sgr A* (jMuno et al.ll2005 ): these 
appear to be LMXBs produced by three-body interac¬ 
tions between stellar binaries and either black holes or 
neutron stars located in the c entr al p arsec (similar to the 
case in globul ar c l uster cores: iMuno et alJ[2005T ) . Indeed, 
IFaucher-Giguere & Loebi (1201 ill argue that the GC envi¬ 
ronment is likely to produce pulsar-black hole binary sys¬ 
tems, via three-body interactions between stellar black 
holes and recycled pulsar binaries; the resulting pulsars 
would have peri ods in the MSP range . Note that, as 
emphasized bv lFaucher-Giguere fe Loebi (120111) . their es¬ 
timate of ten times more MSPs in the GC region than in 
Terzan5 is likely to only be good to within an order of 
magnitude. Overall, though, it appears quite plausible 
that the GC pulsar population is dominated by recycled 
MSPs. 

A possible caveat to the above argument is that glob¬ 
ular clusters are not undergoing active star formation. 
As a result, neutron stars in globular clusters are “old” 
systems; non-recycled pulsars would have long-since spun 
down, so the only remaining pulsars are likely to be recy¬ 
cled MSPs. This is an important factor in the dominance 
of MSPs in the globular cluster population. In the case 
of the GC, there is certainly active star formation in the 
vicinity of Sgr A*, and hence there may well be a sizeable 
population of young pulsars. However, the population of 
neutron stars left over from the many earlier generations 
of star formation at the GC is likely to outnumber the 
population from the current generation of star forma¬ 
tion. This earlier population is likely to have been spun 
up, and would now contribute to the present MSP pop¬ 
ulation. In addition, the much higher stellar density at 
the GC than in globular clusters (by about two orders 
of magnitude) implies that the scattering timescales are 
far shorter in the GC environment. Thus, even if there 
is a population of young pulsars in the GC environment, 
they would be spun up to MSPs much more quickly than 
in globular clusters; the result is that MSPs should still 
dominate the GC pulsar population. 

Fig. [T] overlays the sensitivity of the best pres ent 
searches (|Johnston et al.ll2006t iMacauart et all[20101) on 
the luminosity-period distribution of the known pulsar 
population, assuming the “weak scattering” scenario, i.e. 
that GC pulsars are subject to the same temporal smear- 
ing as the GC magn etar. While the 5 GHz search of 
iJohnston et al.l (120061 ) is m ore sensitive to normal p ulsars 
than the 15 GHz search of Macquart et al. (2010), even 
the relatively-benign assumed temporal smearing causes 
a significant reduction in its sensitivity to fast pulsars 
w ith periods < 10 ms. F urther, even the 15 GHz search 
of lMacauart et ahl (12010 1 is sensitive to only a small frac- 
tion (< 4%) of a GC MSP population whose properties 
resemble those of known MSPs. We emphasize that this 
assumes the optimistic scenario in which the temporal 


smearing towards GC pulsars is similar to that measured 
towards the GC magnetar. The dashed curves in Fig. [T] 
show that the situation is even worse if the GC mag¬ 
netar is being seen through a “hole” in the screen, and 
the temporal smearing towards Sgr A* is similar to the 
earlier estimates: the 5 GHz and 15 GHz searches would 
then be entirely insensitive to pulsars of periods < 40 ms. 
Thus, given that the GC pulsar population is likely to be 
dominated by MSPs and that present searches have been 
insensitive to such a population, it appears that an MSP 
population is a viable way of hiding pulsars at the GC, 
and solving the missing pulsar problem. 

In passing, we note that the 2 brightest MSPs in 
Terza n5 would have been de tected in the 15 GHz sur¬ 
vey of lMacauart et all (|2010f ) . if located at the distance 
of the GC, and if the weak scattering scenario indeed ap¬ 
plies to the GC pulsar population. The lack of detections 
in this survey appears surprising if the GC region indeed 
has ~ 1 0 ti mes the num ber of MSPs seen in Terzan5 
(iFaucher-Giguere fe Loebll201ir) . Indeed, this might be 
considered evidence that either the GC sightline is sub¬ 
ject to strong scattering or that the number of MSPs at 
the GC has been over-estimated. However, we caution 
that small number statistics makes it difficult to extrap¬ 
olate from two bright MSPs in Terzan5 to the full GC 
population. It is hence critical to increase the sensitivity 
of GC searches so as to be able to detect a significant 
fraction of the Terzan5 (or known MSP) population, if 
placed at the GC distance. If MSPs remain undetected 
in such searches, the non-detections would essentially im¬ 
ply that either the weak scattering case is not applicable 
to the GC or that the GC pulsar population is somehow 
very different from that in globular cluster environments. 

3. THE DETECTABILITY OF PULSARS AT THE 
GALACTIC CENTER 

In this section, we compute in detail the expected 
signal-to-noise ratio for MSPs at the Galactic Center 
with the most sensitive current and forthcoming radio 
telescopes. These computations are motivated both by 
the arguments of the preceding section concerning MSPs 
at the Galactic Center and by the dramatic increase 
in the_ f racti onal bandwidth of modern pulsar backends 
(e.g. lSiemion et all2013t ). due to which the freq uency de¬ 
pendence of the relevant qu antities (e.g. lLazio fc Cordesl 
119981 IMacauart et al .1120 1 dl) within each observing band 
plays a critical role in determining the optimal observ¬ 
ing frequency. The following analysis takes into account 
variations in the signal-to-noise ratio (S/N) across the 
observing band due to the pulsar radio spectrum, the fre¬ 
quency dependence of the system temperature, and the 
change in the pulse width at different frequencies caused 
by temporal smearing. 


3.1. The S/N for wide bandwidth pulsar searches at the 
Galactic Center 

For a narrow observing bandwidth, Ajz, the S/N for 
the detection of a pulsar of average f lux density S„, pulse 
period P and width W is (jLorimer et al.ll2006f ) 


S/N = 5, 


Gy/n p AvAt Ip — W 


T m 
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where At is the telescope integration time, n p is the num¬ 
ber of polarizations observed, G is the telescope gain 
and T sys is the total system temperature. Note that the 
above equation assumes that S v , T sys and W do not de¬ 
pend on the observing frequency. However, when the ob¬ 
serving bandwidth is large, this assumption breaks down 
and Equation ([1]) must be generalised to incorporate the 
frequency dependence of both the noise and the signal 
across the observing band. The resulting net S/N can be 
estimated by comparing the total signal in the observ¬ 
ing band to the total noise received during the duration 
of the observations. The total signal measured over the 
frequency interval {u\, v 2 ) in a time period At is then 


where T 0 {y) = 350(i//2.7 GHz) -2 7 K dReich et al.1 
I1990H . We further estimate Oq = 0.33° from the 2.7 GHz 
image of lReich et all (|1990T ). We emphasize that this ap¬ 
proximation only applies in the neighbourhood of Sgr A* 
but is sufficient for our purposes, given the high frequen¬ 
cies (> 5 GHz), and hence, the relatively small telescope 
beamwidths under consideration here. The effective con¬ 
tribution of the sky temperature to the system temper¬ 
ature is determined by computing the weighted average 
of the sky emission over the telescope beam, B(0): 


Tg cW 


jT{d)B{G)d 2 G 
f B{6)d 2 6 



rv 2 

S = n p At / S v {v')dv'. (2) 

Next, the noise power per polarization in a single time 
and frequency channel of respective widths St and Sv is 
n = T sys (y)\/ St Sis/G. For observations over a number 
A c h of such time, frequency and polarization channels, 
the noise adds in quadrature to give 


N = 



(3) 


For an observation with n p polarizations, covering a 
frequency interval (i/±, V 2 ), and of duration At, the above 
can be generalized to 


N = 



G 2 (z/) 


W{v') 

P - W{v') 


dv'. 


(4) 


where, in the second equality, we assume a gaussian 
beam shape of full-width-at-half-maximum 6/, = 1.22A/d, 
where d is the telescope diameter or, in the case of an 
interferometer, the diameter of each element of the array. 

We have computed the S/N values for the two most 
sensitive present telescopes covering the frequency range 
~ 5 — 50 GHz, the Green Bank Telescope (GBT) and 
the Karl G. Jansky Very Large Array (VLA). For these 
telescopes, the system temperature (not including the 
contribution from the GC region) at the low GC ele¬ 
vations, the telescope gain and the usable bandwidth 
in different observing bands have been obtained from 
the telescope webpagefl In addition, we also consider 
the proposed sensitivities of the mid-frequency array of 
the Square Kilometre Array Phase-1 (SKAl -MI D) and 
the full Square Kilometre Array (SKA) (|Dewdnev et al.1 

l20H . 

We consider two scattering scenarios: “weak scatter¬ 
ing”, where the pulse smearing timescale is assumed to 
be the same as that measured towards the GC magnetar 
(jSpitler et al.ll2014l ) 


In the above equation, the factor W/(P — W) inside 
the integral represents the fraction of power present in 
the pulsed signal per unit time. 

Next, the telescope system temperature T sys contains 
contributions from the sky and the telescope receiver 
temperature T rec , each of which has different depen¬ 
dences on frequency. The sky temperature itself con¬ 
tains two contributions, one from the bright GC region, 
which we label Tq c> and the other from the atmosphere, 
T a tm- The atmospheric contribution depends on the tele¬ 
scope site and increases with frequency at frequencies 
> 20 GHz. It is convenient to absorb the two telescope- 
dependent contributions, from the telescope electronics 
and the atmosphere above the site, into a single term 
which we label Tr, so that the total system temperature 
is expressed in the form 


T sys (v)=T GG (v)+T R (u) (5) 

The determination of the GC contribution to the sky 
temperature is complicated by the fact that the emission 
near Sgr A* has a strong spatial dependence, so that 
the exact amplitude and spectrum of the contribution 
depends on the telescope beam size (which itself depends 
on frequency). For simplicity, we model the decline of sky 
temperature with angular distance 9 from Sgr A* using 
a Gaussian profile, 


T(8, v) = T 0 {v) exp 


0 2 
A) J 


( 6 ) 


T smear — 1-3 X S > ($) 

and “strong scatteri ng”, where the scatte ring occurs 
130 pc from the GC (|Lazio fc Cordesl fl998T) and, using 
the GC distan ce and angu l ar bro adening size of Sgr 
A* adopted by iBower et all (12014T) . the pulse smearing 
timescale is given by 

Tsmear « 2.1 X 10 2 X S. (9) 

The observed pulse width is the sum in quadrature of 
the pulse scattering timescale, the intrinsic pulse width 
and the detector channel width: 

W = V ( w 50 P) 2 + r 2 meal . + St 2 , (10) 

where the intrinsic pulse width is taken to be a fixed 
fraction, W 50 , of the pulse period. The time resolution of 
the detector, St, is an important factor when the detector 
time resolution is coarser than the intrinsic pulse width. 
In this regime, we see from Equation (fTUl) that, in the 
limit of narrow pulse widths, the observed pulse width 
tends to the detector time resolution, St, and that no 
further improvement in S/N is possible with a further 
decrease in w^qP. We will assume a time resolution of 
50 /xs for the VLA, SKAl-MID and full SKA, which is 

4 https://science.nrao.edu/facilities/vla 

https://science.nrao.edu/facilities/gbt 
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more than sufficient to not hinder the detection of MSPs. 
However, in the case of the GBT, we have assumed a time 
resolution of 0.5 ms, the best available with the present 
wideband VEGAS spectrometer; this does influence the 
detectability of the fastest MSPs, with periods of ~ 1 ms. 

We note that Equation (JTJ, and its generalization for 
large bandwidths, implies that a pulsar whose pulse 
width W matches the pulse period P would be com¬ 
pletely undetectable, with the noise term increasing 
rapidly as W approaches P. When including the effect of 
scatter broadening on the pulse width, Equation (fTOl) im¬ 
plicitly assumes that the scatter-broadened pulse profile 
is smeared into a boxcar of width W, rather than be¬ 
ing convolved with an exponential tail. The latter case 
would result in so me pulsed emission be ing visible even 
with W > P (e.g. iMacouart etahil201 Q) . A more accu¬ 
rate treatment would slightly improve the detectability of 
pulsars at lower frequencies, but would not significantly 
alter our results. 

3.2. Application to the Galactic Center 

Incorporating all the above considerations, the signal- 
to-noise ratio for a pulsar at the Galactic Center is 


S/N 


1.49 x tQ L 1A (t^sP 


1 r «2 

W( v>) 

y Jv 1 G' z (v') 

P-W(v') 



where we have defined the 1.4 GHz pulsar (pseudo- 
)luminosity L 1.4 = S u d,Q C in units of rnJykpc 2 , the 
numerical prefactor assumes a GC distance of dec = 
8.2 kpc, and a is the pulsar spectral index (£„ oc v~ a ). 

Analyses of MSP spectra indi cate a mean spec¬ 
tral index in the range 1.6-1.8 (|Kramer et all fl998l 
[Maron et al.1 l2000h : we hence adopt a typical value of 
a = 1.7 in our modelling. However, we note that a 
recent analysis of the slow pu lsar pop ulation favours a 
lower spectral index, a = 1.4 dBates et al.ll2013 ). If this 
also applies to the MSP population, it would assist the 
detectability of pulsars at high frequencies. 

The optimal detection frequency for a pulsar of a given 
period is chiefly determined by the competition between 
scattering, which broadens the pulses to a timescale less 
than the pulse period only above some frequency t'Hm, 
and the pulsar spectrum, which declines with increasing 
frequency (of course, secondary considerations include 
the sky contribution to the system temperature, the tele¬ 
scope characteristics and the observing bandwidth). It 
is obvious that the pulsar is undetectable at frequen¬ 
cies below which the pulse broadening timescale matches 
the pulsar period. It might therefore be supposed that 
the optimal detection frequency includes frequencies just 
above uu m . However, the noise contribution becomes ar¬ 
bitrarily large for frequencies in the vicinity of v u m , where 
P mW. In fact, there is a range of frequencies above ou m 
where W[y) is not sufficiently small compared to P, and 
where the contribution of the noise power per frequency 
interval in Equation (Hill exceeds the contribution from 
the pulsed signal power, so that inclusion of this region 
decreases the S/N. It is therefore advantageous to restrict 
the lower bound of the observing band to exclude this re¬ 
gion, so as to maximize the S/N. In our calculations of 


the S/N for bands near ^n m , we therefore optimize the 
S/N by restricting the lower cutoff frequency to include 
only those frequencies which make a positive contribu¬ 
tion to the S/N (i.e. we do not necessarily use the entire 
available bandwidth in the detection). Observationally, 
one would similarly attempt to maximize the S/N of any 
detection by restricting the range of frequencies only to 
those in which pulsed power is evident. 

The four panels of Figure [2] illustrate the detectabil¬ 
ity of MSPs of three different periods, 1 ms, 5 ms and 
20 ms (all with w$o = 10% and Pi .4 = 10 rnJykpc 2 , 
and located at the Galactic Center distance), with a 30- 
hour integration with the four telescopes under consider¬ 
ation. The filled squares and solid curves correspond to 
the low-scattering magnetar-like scenario (with r smea r = 
l-S^GHzS)) while the open circles and dotted curves cor¬ 
respond to the high-scattering scenario arising from a 
screen located 130 pc from SgrA* (T smea r = 208 ^q^ z s). 

It is clear from the figure that, in the weak-scattering 
case, the peak S/N is obtained at a central frequency of 
ss 10 GHz (i.e. the X-band for the GBT and the VLA) 
for MSPs with periods < 10 ms, while for slower (20 ms) 
pulsars, the peak S/N shifts to ~ 6 GHz (i.e. the C-band 
of the GBT and the VLA), and to progressively lower 
frequencies for progressively slower rotators. Conversely, 
in the high-scattering case, the peak S/N for MSPs is 
obtained at much higher frequencies, ss 20 — 30 GHz for 
MSPs with periods in the range 1 — 20 ms. 

In passing, we note that the large data volumes im¬ 
ply that a high statistical significance (typically, > 10 cr) 
is usually required in such pulsar surveys. This implies 
that an MSP with Pi .4 = 10 rnJykpc 2 and with periods 
< 20 ms would be detectable with integration times of 
« 10 — 30 hours with the GBT and the VLA (and, of 
course, with the SKA1-MID and the full SKA) in the 
weak-scattering case. However, much larger integration 
times (> 100 hours) would be needed on the GBT, the 
VLA or the SKAl-MID to detect such an MSP in the 
high-scattering case. If the GC pulsar population is in¬ 
deed dominated by MSPs and the high-scattering case 
applies, it may only be possible to detect this population 
with the full SKA. 

Of course, the optimal detection frequency for MSPs 
critically depends on both the MSP period distribution 
and the MSP luminosity function. Assuming that these 
are same for the GC environment as those o f the known 
MSP population (e.g. lManchester et al.ll2005l h we can de¬ 
termine the overall optimal detection frequency for MSPs 
(integrated over all spin periods), by estimating, for each 
spin period, the frequency at which the peak S/N occurs, 
and then plotting this frequency against the peak S/N 
weighted by the relative fraction of MSPs with this pe¬ 
riod. The results of this analysis for the GBT, VLA and 
SKAl-MID are shown in Figure [H in both the weak- and 
strong-scattering regimes. We note that the period dis¬ 
tribution of the known MSP population is dominated by 
spin periods between 1 and 3 ms. In the weak-scattering 
regime, the S/N for the three telescopes is either rela¬ 
tively flat over the frequency range 5 — 10 GHz or peaks 
at « 10 GHz. It is hence not surprising that the optimal 
detection frequency in the weak-scattering regime closely 
matches that at which the S/N is maximal for pulsars of 
those periods (i.e. « 8 — 10 GHz). However, we note 
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Fig. 2.— The frequency dependence of the S/N for GBT, VLA, SKA-MID and full SKA searches for MSPs at the Galactic Center. The 
filled and open symbols are, respectively, for the weak-scattering (i.e. magnetar-like) and strong-scattering scenarios. The symbols plotted 
are for MSPs with the listed periods (1, 5, and 20 ms), a duty cycle of 10% and a spectral index of —1.7. The vertical axis gives the expected 
S/N after a 30-hour integration for an MSP with a 1.4GHz luminosity of L 1.4 = lOmJy kpc 2 . The frequency ranges of the current suite 
of available receivers for each telescope is shown at the top of each panel. The optimal frequency band for MSPs with periods < 10 ms is 
8 — 12 GHz in the weak-scattering case, and 18 — 26 GHz (i.e. K-band) in the strong-scattering case. 


that in the strong-scattering regime, the optimal detec¬ 
tion frequency is actually ~ 15 GHz for the VLA, i.e. 
lower than the frequency at which the S/N is maximal 
for MSPs with periods of 1 — 3 ms. This is because the 
VLA detection S/N in the strong-scattering regime is far 
higher for MSPs with periods of ss 20 ms than for MSPs 
with periods of ~ 1 — 3 ms, and this compensates for the 
larger fraction of fast MSPs in the population. 

We emphasize that the above analysis makes the crit¬ 
ical assumption that both the MSP luminosity function 
and the MSP period distribution in the GC environment 
are the same as those of the known MSP population, 
dominated by MSPs in globular clusters. It is not im¬ 
plausible that these are different in the GC and globular 
cluster environments. Interactions that form MSPs in 
the GC region are likely to differ substantially from those 
that operate in other known environments (e.g. due to 
the inc reased p revalence of b inar y inte ractions in the GC 
region; iFaucher-Gieuere fe Loebll201lD . 

Recognising that X-band is the optimal detection band 
for most GC MSPs in the weak-scattering regime, Fig¬ 


ure [4] shows the limiting detection luminosity, scaled 
to 1.4 GHz, for a 30-hour survey of the Galactic Cen¬ 
ter as a function of pulse period. The 1.4 GHz lumi¬ 
nosities of_^nowu__pulsars (from the ATNF pulsar cat¬ 
alogue; iMancfresteLjEld] [2005;) are plotted for compar¬ 
ison. As noted earlier, previous pulsar surveys of the 
GC have been insensitive to all but the most luminous 
members of the known MSP population, if located at 
the GC distance: th e most sensitive previous search 
l|Macauart et all f2010H would have detected only 4% of 
all known pulsars with periods < 20 ms. By contrast, a 
30 h X-band survey with the GBT would detect 27%, and 
a corresponding VLA survey 42%, of this population. Of 
course, both such surveys would, in addition, be sensitive 
to > 80% of slow pulsars, with periods > 0.1s. 

3.3. Searches for accelerated pulsars 

Implicit in the foregoing calculations and time esti¬ 
mates is the assumption that it is possible to integrate 
coherently for several tens of hours on individual pul¬ 
sars. This is potentially an issue for MSPs, which tend 
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Fig. 3. — An estimate of the optimal detection frequency for 
MSPs with spin periods up to 50 ms, under the assumption that 
the MSP period distribution at the GC is the same as that of the 
known MSP population. For each frequency band, the plot shows 
the fraction of pulsars with spin periods whose peak S/N falls in 
this band, weighted by their peak detection S/N. The filled symbols 
designate values for the weak scattering case, while open symbols 
indicate the corresponding values for the strong scattering case. 

to occur in binary systems and where the inherent orbital 
accelerations can lead to a non-negligible pulse phase 
drift for observations lasting longer than a small frac¬ 
tion of the orbital period. For example, more than half 
of the known MSPs (with periods < 20 ms) in globu¬ 
lar clusters arise in binarie^l, with orbital peri ods r ang¬ 
ing from » 1.6 hours to ~ 191 days (e.g. ICamilo et al.1 
120001 : iRansom et al.ll2005lh and a median orbital period 
of ss 17 hours. Searches for MSPs at the Galactic Center 
that extend for more than a few hours are hence likely 
to be affected by the above issue, which is relevant to 
about half the MSP population (assuming that the dis¬ 
tribution between “isolated” and “binary” MSPs is the 
same as that in globular clusters). Pulse phase drift is 
a particularly important consideration for GC searches 
with the GBT and VLA, whose northern latitudes re¬ 
strict them to individual GC observing runs of durations 
< 6 hours. A 30-hour GBT or VLA integration on the 
GC would thus require that one integrate coherently over 
a period of at least 5 days. 

Drifts in the pulse period over the duration of the ob¬ 
servation must be taken into account, so that the sen¬ 
sitivity is not reduced for MSPs in short-period bina¬ 
ries. When determining the optimal method for taking 
the pulsar orbit into account in the detection process, 
it is convenient to characterise binary pulsars in terms 
of the ratio of total observation duration, AT, to the 
orbital period, T or b- For pulsars with AT/T or b < 0.1, 
the change in pulse frequency is sufficientl y small that 
standard acceleration sear ches are efficient dWood et al.l 
I199H IRansom et al.1 [20021') . Conversely, if AT/T orb > 
1.5, phase modulation searches may be empl oyed at the 
pena l ty of a slight decreas e in sensitivity (|Jouteux et al.1 
120021 IRansom et al.l I2003T ). However, in the intermedi¬ 
ate regime 0.1 < AT/T or b < 1.5, full sensitivity searches 
require a search over the Keplerian parameters of the or¬ 
bit. For orbital periods similar to those typical of binary 
MSPs in globular clusters, searches over the orbital pa¬ 
rameters would be needed for integrations lasting longer 

5 See http://www.naic.edu/~pfreire/GCpsr.html 


than ss 1 hour. 

In the case of SKAl-MID and the full SKA, the rel¬ 
atively high sensitivity implies that short searches (of 

< 1 hour duration), in conjunction with existing acceler¬ 
ation search software, would be capable of detecting sig¬ 
nificant fractions of the known MSP population (a: 25% 
for SKAl-MID and ~ 90% for the full SKA) if located at 
the GC distance. However, for the GBT and the VLA, 
such short integrations would only be able to detect a 
small fraction (a few %) of the population, as the lim¬ 
iting luminosity would be a factor of « 5.5 worse than 
that of the 30-hour integrations shown in Fig. 01 One 
might increase the net statistical significance of detec¬ 
tions by performing multiple searches over multiple short 
stretches of data. In such an approach, one would re¬ 
tain all candidates detected at a lower significance (i.e. 

< 10cr) and then jointly analyse all such marginal candi¬ 
dates to test whether the different observing epochs yield 
pulse periods and accelerations consistent with a fitted 
orbit. 

The alternative approach would be to fit different com¬ 
binations of the six Keplerian orbital parameters to the 
phase-coherent pulse data. The most effective current 
search methods in this regime have been developed and 
successfully i mplemented to detect binary pulsar s on cir¬ 
cular orbits dKnispell 120111 : iKnispel et al.l I2013T) . This 
technique relies on convolution of the data with multiple 
orbital parameter combinations; its extreme computa¬ 
tional expense has necessitated implementation on highly 
distributed machines (e.g. Einstein@Home). A potential 
disadvantage of this approach is the increase in the di¬ 
mensionality of the search over a conventional acceler¬ 
ation search, requiring a commensurate increase in the 
S/N of any detection in order to achieve an acceptable 
level of significance. However, this can be partially miti¬ 
gated by the restriction of the search to a specific range 
of orbital periods and to small eccentricities, as has been 
done in previous searches. 

4. CONCLUSIONS 

The high stellar density in the central parsec of the 
Milky Way offers a solution to the long-standing mystery 
of the lack of pulsar detections at the Galactic Center. 
Analogous to the situation in globular clusters (where the 
densities are a few orders of magnitude lower than at the 
GC), the high density is likely to result in a GC pulsar 
population dominated by millisecond pulsars. Previous 
GC pulsar searches have been almost entirely insensitive 
to such an MSP population. 

We estimate the optimal search frequency for an MSP 
population for two present and two future telescopes, the 
VLA, the GBT, the SKAl-Mid and the full SKA, as¬ 
suming that the GC pulsar population has a luminosity 
distribution similar to that of field pulsars. We consider 
two scattering cases, weak-scattering, where the scatter¬ 
ing screen is roughly midway between us and the GC, 
and strong-scattering, where the screen is « 130 pc from 
the GC. We find that the optimal MSP search frequencies 
are ~ 8 GHz and ~ 25 GHz in the weak-scattering and 
strong-scattering cases, respectively. Deep (10 —30 hour) 
integrations with the VLA or the GBT should allow a de¬ 
tection of MSPs at the GC at > 10cr significance if the 
weak-scattering case is indeed applicable, as suggested by 
the recent detection of a magnetar close to the GC. How- 





























ever, the strong-scattering scenario would require the full 
SKA to detect and time MSPs at the distance of the 
Galactic Center. 
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Fig. 4. — The 10cr sensitivities of GBT, VLA, SKA1-MID and 
SKA 30-hour X-band integrations to the known Galactic pul¬ 
sar population, if placed at the distance of the Galactic Cen¬ 
ter, assuming the weak-scattering case. As in Fig. [I] the dots 
show the 1.4 GHz lum inosity of the known pulsar population 
([Manchester et ah! I2005T) plotted versus pulsar period, while the 
solid, dashed, dotted and dash-dotted curves show the 10<r sensi¬ 
tivities for the GBT, VLA, SKA1-MID and full SKA, respectively. 
It is clear that deep X-band observations with existing telescopes 
(the GBT and the VLA) would be sensitive to a significant fraction 
(> 30%) of the known MSP population (as well as to > 65% of the 
entire known pulsar population), if located at the GC distance. 
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